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Hysteretic properties of a magnetic particle with strong surface anisotropy 
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We study the influence of surface anisotropy on the zero-temperature hysteretic properties of a 
small single-domain ferromagnetic particle, and investigate limiting cases where deviations from the 
Stoner-Wohlfarth model are observed due to non-uniform reversal of the particle's magnetization. 
We consider a spherical particle with simple cubic crystal structure, a uniaxial anisotropy for core 
spins and radial anisotropy on the surface. The hysteresis loop is obtained by solving the local (cou- 
pled) Landau-Lifshitz equations for classical spin vectors. We find that when the surface anisotropy 
constant K s assumes large values, e.g. of the order of the exchange coupling J, large deviations are 
observed with respect to the Stoner-Wohlfarth model in the hysteresis loop and thereby the limit- 
of-metastability curve, since in this case the magnetization reverses its direction in a non-uniform 
manner via a progressive switching of spin clusters. This characteristic value of K 3 depends on the 
surface-to-volume ratio of exchange coupling and the angle between the applied field and core easy 
axis. 

PACS number(s): 75.50.Tt - 75.30.Pd - 75.10.Hk 
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I. INTRODUCTION 

Surface effects have a strong bearing on the proper- 
ties of small magnetic systems, and entail large deviations 
from the bulk behavior. It was shown in 1 that the mag- 
netic disorder on the surface caused by surface anisotropy 
is long ranged, which implies that even the spins in the 
core of a very small magnetic particle (2 nm) render a 
magnetization that deviates from the bulk value. It will 
be useful to understand surface effects in magnetic mate- 
rials in order to control their properties which are relevant 
for technological applications. One such property is the 
coercive field as it gives indications on the relaxation time 
of the magnetization and thereby on the stability of the 
information stored on magnetic media. 

Surface effects are due to the breaking of crystal- 
field symmetry, and this is a local effect. So, in order 
to study such effects one has to resort to microscopic 
theories, unlike the macroscopic Stoner-Wohlfarth (SW) 
model 2 , which are capable of distinguishing between dif- 
ferent atomic environments and taking account of phys- 
ical parameters such as bulk and surface anisotropy, ex- 
change and dipole-dipole interactions. Unfortunately, 
this leads to difficult many-body problems which can only 
be dealt with using numerical approaches. 

This work deals with the effect of strong surface 
anisotropy on the hysteretic properties (hysteresis loop 
and limit-of-metastability curve, the so-called SW as- 
troid), of a single-domain spherical particle (with free 
surfaces), a simple cubic (sc) crystal structure, a uniaxial 
anisotropy in the core, and radial single-site anisotropy 
for spins on the boundary. The hysteresis loop and 
thereby the critical field are computed by solving, at zero 



temperature, the local Landau-Lifshitz equations derived 
from the classical anisotropic Dirac-Heisenberg model in 
field, subjected to a local condition (see below) account- 
ing for the minimization of energy with respect to local 
rotations of each spin in the particle. In Ref. 3 the same 
method was used for studying the hysteretic properties 
of models of nanoparticles, where the anisotropy was ei- 
ther random in the whole particle or taken only on the 
surface, and the analysis was restricted to the hysteresis 
loop. 

In this paper, we use an improved version of the 
method mentioned above including a global-rotation con- 
dition on the resultant magnetic moment of the particle 
in addition to the local condition (see 4 ). We compute the 
hysteresis loop and infer from it the limit-of-metastability 
curve (SW astroid), and compare with the SW model es- 
pecially when the surface anisotropy constant assumes 
large values, e.g. K s /J ~ 1. This study has allowed us 
to investigate the limit of validity of the SW model for 
very small magnetic particles where surface anisotropy 
plays a determinant role, and whose magnetization no 
longer switches in a coherent way. 

Our method, based on the numerical solution of the 
Landau-Lifshitz equation at zero temperature, is checked 
against the SW semi-analytical results in two limiting 
cases of the exchange coupling with different distributions 
of anisotropy axes. We first consider a single-domain par- 
ticle with a macroscopic magnetic moment resulting from 
very strong exchange interaction. This is equivalent to 
the SW one-spin problem with uniaxial anisotropy. A 
second test deals with the case of a square particle of 
non-interacting spins all with randomly distributed easy 
axes. This model mimics an assembly of mono-dispersed 
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single-domain nanoparticles with a random distribution 
of their easy axes embedded in a 2d non-magnetic matrix. 

The plan of this work is as follows: we first define our 
model (Hamiltonian and physical parameters), present 
the method used for computing the hysteresis loop, and 
test it against the semi-analytical results of SW model. 
Then, we discuss our results for a spherical particle in 
terms of exchange coupling, particle's size, and surface 
anisotropy by varying, in turn, one of them while keeping 
the other two fixed. We also study the situation with (in- 
tra) surface exchange coupling different from that in the 
core of the particle. A short account of the present work 
can be found in Ref. 5 . It is worth mentioning though that 
in fact only anisotropy and exchange coupling on the sur- 
face can be considered as free parameters as there are so 
far no definite experimental estimations thereof. 



II. MODEL HAMILTONIAN 

We consider the following classical anisotropic Dirac- 
Heisenberg model 



(1) 



where Sj is the unit spin vector on site i, H is the uniform 
magnetic field applied in a direction ip with respect to the 
reference z axis, Af is the total number of spins (core and 
surface), and in the sequel D will denote the particle's 
diameter. J t j (— J > 0) is the strength of the nearest- 
neighbor exchange interaction, which will be taken in our 
calculations the same everywhere inside the particle, un- 
less otherwise specified (see Fig. 14 et seq); H an is the 
uniaxial anisotropy energy 
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with easy axis ej and constant Ki > 0. This anisotropy 
term contains either of the two contributions stemming 
from the core and surface, and depends on the system un- 
der consideration. For instance, for the 2d model (which 
serves as a test of our calculations by comparison with the 
SW model) all spins (core and surface) have the same 
anisotropy constant but randomly distributed axes. In 
the case of a spherical particle, all core spins are at- 
tributed the same constant K c and all surface spins are 
attributed the constant K s . Moreover, core spins will 
have an easy axis along the z axis, whereas a surface spin 
is assumed to have its anisotropy axis along the radial 
direction, see 6 and many references therein. 

A more physically appealing microscopic model of 
surface anisotropy was provided by Neel 7 , 



H. 



Neel 



(3) 



where Zi is the coordination number of site i and e^- = 
Yij I rij is the unit vector connecting the site i to its near- 
est neighbors. This model is more realistic since the 



anisotropy at a given site occurs only when the latter 
looses some of its neighbors, i.e. when it is located on the 
boundary. However, the extra sum on nearest-neighbors 
in (3) makes this model less practical for numerical cal- 
culations, especially those that are time consuming, such 
as the SW astroid. So in this paper we restrict ourselves 
to the model of radial single-site anisotropy on the sur- 
face. In Ref. 8 , we have developed an analytical theory, 
together with the numerical method used here, for weak 
surface anisotropy and studied this model and compared 
it with the radial-anisotropy model. 

A remark is in order concerning the dipole-dipole in- 
teractions inside the particle. It is well known 9 that these 
relativistic interactions lead to two contributions, a first 
term that is an integral over the volume of the particle, 
and a second one over the surface. The latter represents 
the magnetostatic energy. However, it has been shown 10 
that in very small particles the first contribution is negli- 
gible as compared with the contribution of exchange in- 
teractions. On the other hand, the second contribution 
plays the role of shape anisotropy, which for a spherical 
particle yields an irrelevant constant. Therefore, in our 
case of very small spherical particles, where the effect of 
surface anisotropy constant is most important, which is 
one of the main issues of the present work, the volume 
term is negligible and the shape anisotropy is absent. 



III. METHOD OF CALCULATION OF THE 
HYSTERESIS LOOP 

Different models of a nanoparticlc are studied. In 
each case, we simulate the lattice with sc crystal struc- 
ture, and then assign to each site a length-fixed three- 
component spin vector. For the calculation of the hys- 
teresis loop we start with a magnetic configuration where 
all spins are pointing in the same direction —z, which 
corresponds to the saturation state. The hysteresis loop 
is due to the existence of metastable states in the system. 
Starting from the initial configuration and applied field, 
the integration of the Landau-Lifshitz equation (see be- 
low) tends towards a new configuration that is an energy 
minimum. 

Let us now establish the Landau-Lifshitz equations 
for the magnetic moments. We choose K c as the energy 
scale and normalize the other physical constants accord- 
ingly, i.e., 
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Then, the Landau-Lifshitz (LL) equation for a spin Sj at 
site i, reads 



J = -S t x hf f - aS t x (St x h! 



eff 



(5) 



where a(~ I) is the damping parameter and is the 
effective field acting on the spin Sj and is given by 
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h f / = h +^X>A+ h r (6) 

3 = 1 

where hf" = -{dH an /dSi)/2K c , with ff a „ given in Eq. 
(2), Zi is the coordination number of site i. In the se- 
quel, we will use the reduced parameters, j = J/K c , 
k s = K s /K c . Therefore, for each site i we arrive at three 
coupled equations (for Sf, Sf , Sf), and because of the sec- 
ond term in (6) we actually obtain a system of 37V (local) 
coupled equations. We emphasize that it is more conve- 
nient to use spherical coordinates (for each spin) instead 
of the Cartesian ones. Indeed, owing to the fact that 
the spins are of constant length, this reduces the number 
of individual (for each spin) equations to two instead of 
three, 

9 i = {h%"+ah%") i (7) 
<Pi = (-htff + ah'K'jjsmei, 

where h e J^ = —dH/d9,h e Jl = —dH/dip, are the po- 
lar components of the effective field. For the one-spin 
problem, these are obtained by direct differentiation of 
the energy written in spherical coordinates, whereas for 
a particle it is not possible to obtain a tractable analyt- 
ical expression of the energy in spherical coordinates, so 
hg^ and h e J^ are written in terms of the time derivatives 

of the Cartesian components of in (6). Using Eq. 
(7) instead of (5) allows for a gain of computer time, but 
this method encounters stability problems specific to the 
spherical coordinates, because of the factor 1/ sin0 in (7), 
which diverges as 6 — > 0,ir, and hence a special care is 
required when numerically handling these equations. 

After having constructed the magnetic structure (lat- 
tice and spin vectors on it), we apply a magnetic field H 
at some angle ip with respect to the reference z axis, with 
values chosen in a regular mesh. Then we calculate the lo- 
cal effective field for all spins and thereby the right-hand 
sides of the LL equations (7) and proceed with the time 
integration. As this is done, the total energy in Eq.(l) 
smoothly decreases, and some criterion must be used for 
stopping the integration for the given value of the applied 
field and moving to the next value. In our calculations 
we proceed to the next field value when 

1 dSi 

which implies that the system is close to a station- 
ary state, e being a small parameter of the order of 
1CP 5 — 10~ 7 . However, it was shown in 4 , that this lo- 
cal condition, which accounts for the minimization of en- 
ergy with respect to local rotations (or small deviations) 
of each spin, must be supplemented by a global condition 
on the resultant magnetic moment so as to account for the 
global rotation of the particle's magnetic moment. Obvi- 
ously, for a single spin these two conditions boil down to 
one and the same condition (8). 



Next, the stationary state thus obtained is used as 
the initial state for the next value of the field. Iteration 
of this process over a sequence of applied fields, of given 
magnitude and direction ip, renders the hysteresis loop. 
For each value of this angle we determine the critical or 
switching field (see discussion below). The whole pro- 
cedure finally renders the critical or switching field as a 
function of the angle ip, which in the case of critical field 
is the SW astroid. 

As a test of this method, we considered a box- 
shaped particle with 11 TV = 3 3 , a sc structure, uniax- 
ial anisotropy, and strong exchange interaction between 
spins inside the particle, and computed the hysteresis 
loop for different values of the angle ip between the ap- 
plied field and the easy axis. The results are shown in 
Fig. 1 (left). Next, we present in Fig. 1 (right) the SW as- 
troid, which separates the region with two minima of en- 
ergy from that with only one minimum. We see that the 
SW results are exactly reproduced by our calculations. 
We have also computed the hysteresis loop of a square 
particle of non-interacting spins (J = 0) all with ran- 
domly distributed easy axes. This is equivalent to an as- 
sembly of mono-dispersed single-domain non-interacting 
particles with randomly distributed easy axes in two di- 
mensions. As expected, we find that the remanent mag- 
netization is equal to 1/2. For later reference, we plot in 
Figs. 2 the critical field h c and the height of the magne- 
tization jump (i.e. m u — mj), as functions of the angle 
ip between the direction of the field and core easy axis. 
Obviously, h c (ip) in Fig. 2 (left) is a well known result 
of the Stoner-Wohlfarth model. On the other hand, we 
note that the height of magnetization jump has an al- 
most linear dependence on ip, except for the final portion 
76° < ip < 90°, which corresponds to cycles with cross- 
ing branches as exhibited by the hysteresis for ip = 85° 
in Fig. 1 (left), see 12 for a discussion of this issue. 

IV. SPHERICAL PARTICLES: RESULTS AND 
DISCUSSION 

Here we consider a single-domain spherical particle of 
simple cubic (sc) structure with uniaxial anisotropy in the 
core and anisotropy constant K c , and radial anisotropy 
on the surface with constant K s . Our main goal here 
is to investigate the influence of surface anisotropy, both 
in direction and strength, on the hysteresis loop and SW 
astroid. However, we will also study the effect of exchange 
coupling and particle's size. Again for later reference, we 
plot in Fig. 3 the distribution of surface anisotropy axes 
of the spherical particle as a function of the azimuthal 
angle tp s between a surface spin easy axis and applied 
field. 



A. Effect of the exchange coupling j 

Now we study the effect of exchange coupling on the 
hysteresis loop of a spherical particle containing Af = 360 
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FIG. 1. Left: (numerical) hysteresis loops for different values of ip increasing inwards: ip = 0, 60°, 85°, 90°, for a 3 particle 
with uniaxial anisotropy. For the sake of clarity the SW analytical hysteresis loops have been omitted, since they exactly 
coincide with the computed ones. Right: (numerical in squares and analytical in full line) SW astroid for the same particle; 
J = 10. h 
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FIG. 2. One-spin problem. Left: critical field as function of ip. Right: height of magnetization jump as function of ip. 



spins (176 surface spins and 184 core spins). We first 
consider the case in which the anisotropy constants in 
the core and on the surface are equal, i.e. k s — 1.0, and 
the magnetic field applied along the easy axis of the core 
spins, so as to investigate the influence of radial direc- 
tion of surface anisotropy. For j <C 1, i.e. j = 0,0.01, 
we can see along portion 1-2 in Fig. 4 a progressive de- 
crease (in absolute value) of the magnetization, which is 
due to the alignment of surface spins, since as the field 
direction is along the core easy axis the core spins have 
a rectangular cycle and the jump is at h = 1.0. Next, 
along portion 2-3 we can see two jumps. Indeed, accord- 
ing to the distribution of surface easy axes in Fig. 3, and 
the critical field as a function of ip in Fig. 2 (left), those 
surface spins with ip s between 0.6 and 1.0 are responsible 



for the first jump, and those with ip s between 0.4 and 
0.6 or 1.0 and 1.2 are responsible for the second jump. 
Next, along portion 3-4 we have successive small jumps 
and thereby a slight decrease of the magnetization. The 
origin of these small jumps resides in two contributions. 
One contribution comes from those surface spins whose 
easy axis makes an angle around 0.2 with the field. Even 
though the corresponding height of jump is large (see Fig. 
2, right), their number is rather small (see Fig. 3) thus 
rendering a small contribution to the magnetization. The 
other contribution is due to surface spins with an angle 
ip s ~ 1.4, which yield a small contribution owing to the 
fact that the height of the corresponding jump is very 
small (see Fig. 2, ip s > 1.2), even though their number is 
relatively large. On the last portion of the lower branch 
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FIG. 4. Hysteresis loop, i.e. plot of the magnetization pro- 
jection on the field direction as a function of the (reduced) 
field h, for ip = 0, k s = 1 and different values of j. J\f = 360. 



FIG. 3. Distribution of surface anisotropy axes versus the 
azimuthal angle ip s for a spherical particle with D = 10 
(TV = 360: 176 surface spins and 184 core spins). 

of the hysteresis in Fig. 4, we see another big jump, which 
is due to the switching of core spins at the field h c = 1.0. 
At last, there is a slow increase of magnetization due to a 
final adjustment of surface spins along the field direction. 
In the present case, the surface fully switches before the 
core (see Fig. 5). 

For j = 0.1, we see that the surface behavior remains al- 
most the same as in the previous cases, whereas the core 
spins now switch cluster-wise as can be seen in the 4 th 
picture of Fig. 5. Indeed, regarding the exchange field as a 
small perturbation of the applied magnetic field, it is clear 
that the core spins located near the surface are subject to 
an effective field whose direction is slightly deviated from 
their easy axis, i.e. the corresponding angle tp is slightly 
different from zero. Now, in Fig. 2 (left) we can see that 
this little deviation in ip produces an important change 
in the switching field. On the contrary, we find that this 
effect is almost absent in what concerns the jumping field 
of surface spins, as can be seen along portion 2-3 in Fig. 
4 upon comparing the loops for j = 0, 0.01 and j = 0.1. 
Indeed, the surface spins responsible for these jumps have 
their easy axes at an angle 0.6 < tp s < 1.0, and hence the 
change in the corresponding critical field is very small (see 
Fig. 2 left). In Fig. 4 we can also see that for j = 0.1, 
i.e. when the exchange energy becomes comparable with 
anisotropy and Zeeman energy, there are more jumps that 
can be attributed to the switching of different spherical 
shells of spins starting from the surface down to the cen- 
ter. This situation is sketched in Fig. 5. For example, for 
h = one can see that the exchange has a little influence 



on surface spins, as they are directed almost along their 
easy axes; for h = 0.64 the surface spins show the same 
behavior as in the absence of exchange, but part of core 
spins, located near the surface, are deviated from their 
easy axes. At the field h = 0.8 all these core spins have 
already switched. 

For j = 1 ~ k s , even that there is only one jump, 
the hysteresis loop is not rectangular owing to the fact 
that the spins rotate in a non-coherent way, as can be 
seen in Fig. 6. This is due to a compromise between 
anisotropy and exchange energies, see for example the 
picture for h — 0. Moreover, even a small number of 
neighbors lying in the core produces a large effect via 
exchange on the behavior of a surface spin. For much 
larger values of j the spins are tightly coupled and move 
together, and the corresponding (numerically obtained) 
critical field h c coincides with the (analytical) expression 
obtained in the limit J — > oo, i.e. h c — N c /J\f, where 
N c is the number of core spins. This expression for h c 
has been obtained by summing over the direction of sur- 
face easy axes which results in a constant surface energy 
contribution proportional to k s . Hence, due to spherical 
symmetry, the surface anisotropy constant does not enter 
the final expression of h c . 

Now we consider the case of larger values of k s , e.g. 
k s = 10, so as to investigate the influence of surface 
anisotropy both in direction and strength. The results 
are presented in Fig. 7 (left). Here, a notable difference 
with respect to the previous case, k s = 1, is the fact that 
the core now switches before the surface and at higher 
fields. Moreover, there appear more jumps which may be 
attributed to the switching of various clusters of surface 
spins. Both cases show that as the ratio j/k s decreases, 
the magnetization requires higher fields to saturate. This 
is further illustrated by Fig. 7 (right) where k s = 10 2 = j 
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FIG. 5. Magnetic structure for j — 0.1, k 3 — 1 for the field values h = —4.0,0,0.64,0.8,0.88,4 which correspond to the 
saturation states and different switching fields shown in Fig. 4. These field values correspond to the pictures when starting 
from the upper array and moving right, down left, and then right. Obviously, grey arrows represent core spins and black arrows 
represent surface spins. 



for a smaller particle. 

Let us now summarize the ongoing discussion. We 
observe that considering a radial distribution for sur- 
face anisotropy, leads, even in the case of very strong 
exchange, to an important quantitative deviation from 
the classical SW model. In particular, the critical field in 
our model is given by 

k = jfH:, ( 9 ) 

where H r c is the critical field for a spherical particle with 
radial anisotropy on the surface and uniaxial in the core, 
if" is the critical field for a spherical particle with uni- 
axial anisotropy for all spins. Therefore, when j and k s 
are comparable, the compromise between exchange cou- 
pling, favoring a full alignment of the spins along each 
other, and surface anisotropy, which favors the alignment 
of spins along their radial easy axes, produces large de- 
viations from the SW model. More precisely, the shape 
of the hysteresis loop is no longer rectangular and there 
appear multiple jumps. The appearance of these jumps 
makes it necessary to define two field values with the help 
of which a hysteresis loop can be characterized. A value 
that marks the limit of metastability, called the critical 
field, and the other value which marks the magnetization 
reversal, i.e. when the projection of the magnetization 
on the field direction changes sign, and this is why it is 
called the switching field (or still coercive field). 



B. Effect of the particle's size M 

Here, we study the effect of varying the particle's size 
while keeping j and k s fixed. So we use the same value 
of anisotropy constant for all spins and strong exchange, 
i.e. k s = 1, j = 10 2 , and vary the particle's diameter 
from 6 {Af = 56) to 30 (TV =' 12712). In Fig. 8 (left) 
are presented hysteresis cycles of a particle with different 
diameters when the field is along the core easy axis, and 
on the right the variation with the particle's diameter 
of the critical field 13 (in diamonds) obtained from the 
numerical solution of the Landau-Lifshitz equation for 
j = 10 2 , and (in circles) the SW critical field multiplied 
by the core-to- volume ratio (see Eq. (9)). The figure on 
the left shows that for such a value of k s the hysteresis 
loop is rectangular for all sizes, and that the critical field 
decreases with the particle's size. The latter fact is clearly 
illustrated by the plot on the right, which also shows that 
for k s = 1 = 10~ 2 j, all these hysteresis loops can be 
scaled with those rendered by the SW model. Next, Fig. 
9 shows the variation with the surface-to-volume ratio 
N st = N s /N of the critical field for all angles between 
the core easy axis and magnetic field, this is the limit- 
of-metastability curve. These results show that, even in 
the general case of a field applied at an arbitrary angle 
with respect to the core easy axis, the critical field of a 
spherical particle with k s — 1 can be obtained from the 
SW model through a scaling with constant N c /N. One 
should also note that the astroid for all particle sizes falls 
inside that of SW, in accordance with Fig. 8 (right), and 
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FIG. 7. Left: Hysteresis loops for tj) = 0, fe s = 10 and different values of j. D = 10 (TV = 360). Right: Hysteresis loops for 
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the larger the surface contribution the more the astroid 
shrinks. 

Therefore, for k s = 1 our results for the hystere- 
sis loop and limit-of-metastability curve can be scaled 
with those of SW model with the scaling constant N c /Af, 
which is smaller than 1 for a particle of any finite size. 

Next, in Fig. 10 (left) we present the hysteresis loop in 
the case where the surface anisotropy constant k s equals 



the exchange coupling and the field is applied along the 
core easy axis, and in Fig. 10 (right) the switching field 14 
as a function of the particle's diameter D. There are 
two new features in comparison with the previous case 
of k s = 1: the values of the switching field are much 
higher, and more importantly, its behavior as a function 
of the particle's size is opposite to that of the previous 
case. Indeed, here we see that this field increases when 
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FIG. 9. Astroid for k B — l,j — 10 2 for different values of 
the surface-to- volume ratio N st = N a /M. The lines on the 
astroids inside the SW one are only guides for the numerical 
data. 

the particle's size is lowered. For such high values of k s 
(K s 3> K c ) surface spins are aligned along their easy 
axes, and because of strong exchange coupling they also 
drive core spins in their switching process. Thus, the 
smaller is the particle the larger is the surface contri- 
bution, and the larger is the field required for complete 
reversal of the particle's magnetization. This could ex- 
plain the non-saturation of magnetization that has been 
observed in e.g. cobalt particles 15 . 



C. Effect of the surface anisotropy constant k s 

Now, we fix the exchange coupling constant j, the 
particle's total number of spins J\f, and vary the surface 
anisotropy constant k s . Because K c is in general 2 to 3 
orders of magnitude smaller than J, we have investigated 
the effect of surface anisotropy constant in the case of 
j = J/K c = 10 2 . 

In contrast with the case k s = 1 and j — 10 2 — lO 3 
where the hysteresis loop and the limit-of-metastability 
curve scale with the SW ones with the same scaling con- 
stant for all angles between the applied field and core easy 
axis, we find that for 1 < k s < 20 the scaling constant 
depends on the angle ip, as can be seen in Fig. 11. This 
fact explains the deformation of the SW astroid, that is a 
depression in the core easy direction and an enhancement 
in the perpendicular direction. 

For larger values of k s we have computed the hys- 
teresis loop for ij) = 0,Af = 360, j = 10 2 . The results 
are given in Fig. 12. Here, we first note that the shape 
of the hysteresis loop is rather different from that ren- 
dered by the SW model, since for k s = 30, for instance, 
the hysteresis loop is no longer rectangular, even that 
ijj = 0. As explained earlier, this effect is due to the now 
more pronounced non-uniform rotation of surface spins 
and core spins located near the surface, and thereby that 
of the particle's magnetization. This non-uniform switch- 
ing process causes large deviations from the SW model, 
and thereby no scaling with the latter is possible. From 
Fig. 12, we extract and plot in Fig. 13 the switching field 
h c as a function of k s /j, denoted by k s in the sequel. We 
find that h c first slightly decreases for k s ^0.1 and then 
increases, and when k s approaches 1 it jumps to large 
values. As discussed above, for such high values of k s sur- 



8 




FIG. 10. Left: Hysteresis cycle for if) = 0, j = k s = 10 2 , and different values of the particle's diameter D. Right: Switching 
field as a function of N for the same parameters. 



h 



-► \ 

♦ \ 
♦ \ 


D= 10 (360 spins) 
J-1(f 








o k s = 1 

* k s = 5 

♦ k s = 10 
Scaled SW 




1 1 


•\» 
\ ♦ 
\. ♦ 

1 i 1 i 


♦ 

. 





I -I 1 1 1 1 1 1 1 1 1 — , — 1 I 

0.0 0.1 0.2 0.3 0.4 0.5 h 0.6 



X 



FIG. 11. Astroid for j = 10 2 , Af = 360 and different values 
of surface anisotropy constant k s . The full dark line is the 
SW astroid scaled with N c /M, but the dotted line is only a 
guide for the eye. 

face spins are aligned along their easy axes, and because 
of strong exchange coupling they also drive core spins in 
their switching process, which then requires a very strong 
field to be completed. Clearly, this particular value of k s , 



to be denoted by kg (= 1, here) marks the passage from 
a regime where scaling with the SW results is possible 
(cither with a ^-dependent or independent constant) to 
the second regime where this scaling is no longer possible 
because of completely different switching processes. 

Now we present additional data which show that the 
"critical value" k^ introduced above depends on (at least) 
two parameters. These are the surface-to-core ratio of 
exchange coupling j s / j and the angle ijj at which the field 
is applied with respect to the core easy axis. 

Let us first discuss the effect of (intra) surface ex- 
change coupling. In real materials such as maghemite, 
it was argued in 1 on account of Mossbauer spectroscopy 
that j s /j < 1. In Fig. 14 we have plotted the results for 
h c obtained with surface exchange coupling j s = J s /K c 
smaller than j, i.e. core-core and core-surface couplings. 
First, we see that the "critical" value k c s of k s separat- 
ing the two regimes discussed above, decreases with the 
ratio js/j- This is a consequence of the fact that when 
j s / j < 1, surface spins align more easily along their (ra- 
dial) anisotropy axes since now they experience a weaker 
effective field. We also note that the jump becomes 
smoother. Next, if we consider the curve j s /j = 1 to- 
gether with any other curve with k s /j < 1, we see that 
when k s < 1 the switching field is larger for j s < j than 
for j s = j, and the opposite holds when k s > 1. 

To understand this result, let us imagine a particle 
containing (at least) two groups of surface spins, a group 
1 with exchange coupling j s — j and group 2 with j s < j. 
When k s < k s < k sl k l s being the critical value of k s 
for group i, the spins in group 1 are of SW type, while 
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FIG. 12. Hysteresis loop for ip = 0, j = 10 2 , D — 10 and different values of surface anisotropy constant k s . These two sets of 
data can not be presented as one plot because of scaling mismatch. 



those of group 2 are of non-SW type, in the sense that 
they switch in a coherent way or cluster-wise, respec- 
tively. Hence, as demonstrated earlier, the reversal of 
spins in group 2 always requires a larger switching field. 
On the other hand, when k s exceeds the largest exchange 
coupling in the particle, i.e. j, the switching field of the 
whole particle decreases with j s /j. Now the spins of both 
groups are of non-SW type, and their switching operates 
cluster-wise, but obviously the latter requires a higher 
applied field for group 1 than for group 2. 

Next, a similar effect is obtained when the field is 
applied at an arbitrary angle with respect to the core 
easy axis, as is the case for instance in an assembly of 
nanoparticles. Here, we consider the case of ip = 7r/4. We 
find that there appear multiple large jumps at a smaller 
value of k s (~ 0.2), as can be seen in Fig. 15. 

For an order of magnitude estimate of K s and the 
critical (or saturation) field, consider a 4 nm cobalt par- 
ticle of fee crystal structure, for which the lattice spac- 
ing is a = 3.554 A, and there are 4 cobalt atoms per 
unit cell. The (bulk) magneto-crystalline anisotropy is 
K c ~ 3 x 10~ 17 erg/spin or 2.7x 10 6 erg/cm 3 , and the sat- 
uration magnetization is M s ~ 1422 emu/cm 3 . The criti- 
cal field is given by H c = (2K C /M s )h c . For ip = 0, k c s = 1 
and h c — 15, so H c ~ 6T. On the other hand, k L s = 1 



means that the effective exchange field experienced by 
a spin on the surface is of the order of the anisotropy 
field, i.e. zSJ/2 ~ 2K S . Then using J ~ 8mev we get 
K s ~ 5.22 x 10~ 14 erg/spin, or using the area per sur- 
face spin (approximately a 2 /8), K s ~ 5erg/cm 2 . For the 
case of ip — 7r/4, k c s ~ 0.2 and h c ~ 0.3, which leads to 
H c ~ 0.1 T and K s ~ 1.2 x 10~ 14 erg/spin or 1.2 erg/cm 2 . 

V. CONCLUSION 

Our model of a spherical particle with uniaxial 
anisotropy in the core and radial anisotropy on the sur- 
face leads to mainly two pertinent regions for the surface 
anisotropy constant fc s , with k s > 1 (K s > K c ): 

• For small values of this parameter, e.g. k s /j ~ 0.01 
our model renders hysteresis loops and limit-of- 
metastability curves that scale with the SW results 
for all values of the angle ip between the core easy 
axis and the applied field, the scaling constant be- 
ing N c /Af, which is smaller than 1 for a particle of 
any finite size. On the other hand, the critical field, 
which coincides in the present case with the switch- 
ing field, increases with the particle's size and tends 
to the SW critical field in very large systems, and 
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FIG. 13. Switching field versus the surface anisotropy con- 
stant for ip = 0, j = 10 2 , and D = 10. 

thereby the corresponding astroid falls inside the 
SW astroid for all particle sizes. 

For larger values of k s /j, but k s /j < 0.2, we still 
have the same kind of scaling but the correspond- 
ing constant depends on tp. This is reflected by 
a deformation of the limit-of-metastability curve. 
More precisely, the latter is depressed in the core 
easy direction and enhanced in the perpendicular 
direction. However, there is still only one jump in 
the hysteresis loop implying that the magnetization 
reversal can be considered as uniform. 

• For much larger values of k s j j, starting from k s / j ~ 
1, there appear multiple steps in the hysteresis loop 
which may be associated with the switching of spin 
clusters. The appearance of these steps makes the 
calculated hysteresis loops both qualitatively and 
quantitatively different from those of SW model, as 
the magnetization reversal can no longer be con- 
sidered as uniform, and one has then to define two 
characteristic values of the field associated with a 
hysteresis loop: the critical field and the switching 
field. In addition, in the present case, there are two 
more new features: the values of the switching field 
are much higher than in SW model, and more im- 
portantly, its behavior as a function of the particle's 
size is opposite to that of the previous cases. More 
precisely, here we find that this field increases when 
the particle's size is lowered. This is in agreement 
with the experimental observations in nanoparticles 
(see e.g. 16 for cobalt particles). 
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FIG. 14. Switching field versus the surface anisotropy con- 
stant for ip = 0, and different values of surface-to-core ratio of 
exchange couplings; D = 10. 
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FIG. 15. Hysteresis for ip = 45°, j = 10 2 ,D = W(Af = 360) 
and different values of surface anisotropy constant k s . 

Therefore, assuming radial anisotropy on the surface, 
we find that there is a "critical" value (K s /J) c of the 
ratio K s /J beyond which, large deviations are observed 
with respect to the SW model in the hysteresis loop and 
thereby the limit-of-metastability curve, since in this case 
the magnetization reverses its direction in a non-uniform 
manner via a progressive switching of spin clusters. So, 
in order to deal with these new features one has to resort 
to microscopic approaches such as the one used in this 
work. In fact, it is found that the critical value (K s /J) c is 
even smaller for smaller surface-to-core ratios of exchange 
coupling and larger angles between the applied magnetic 
field and the core easy direction, as it is more likely in 
realistic materials. 
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In a subsequent work we apply the present method 
to cubo-octahedral cobalt particles with a diameter of ca. 
3 nm recently studied in 17 (see also 18 for Pt particles). 
These are particles with fee structure and truncated oc- 
tahedrons on the surface, in which the core has a cubic 
anisotropy, and the surface anisotropy easy axes are be- 
lieved to be along edges and facets with different con- 
stants Kf but whose values are uncertain at present. In 
our calculations we vary these parameters and study the 
effect of surface anisotropy on the Stoner-Wohlfarth as- 
troid that has been experimentally measured in 17 , where 
these anisotropy constants have been estimated from 
magnetic measurements. The final outcome of our cal- 
culations should give an estimation of Kf by comparing 
with these experimental results. Another related issue of 
particular interest to us is the fact that these fee parti- 
cles (see 17 for cobalt and 19 for iron) seem to exhibit an 
effective uniaxial anisotropy despite their cubic crystal 
symmetry. This work is in progress. 
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